Electronic structure of Pr^Ce^CuC^ studied via ARPES and LDA+DMFT+E k 
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The electron-doped P^-iCe^CuCU (PCCO) compound in the pseudogap regime (x « 0.15) 
was investigated using angle-resolved photoemission spectroscopy (ARPES) and the generalized 
dynamical mean-field theory (DMFT) with the k-dependent self-energy (LDA+DMFT+Ek) ■ Model 
parameters (hopping integral values and local Coulomb interaction strength) for the effective one- 
band Hubbard model were calculated by the local density approximation (LDA) with numerical 
renormalization group method (NRG) employed as an "impurity solver" in DMFT computations. 
An "external" k-dependent self-energy Ek was used to describe interaction of correlated conducting 
electrons with short-range antiferromagnetic (AFM) pseudogap fluctuations. Both experimental and 
theoretical spectral functions and Fermi surfaces (FS) were obtained and compared demonstrating 
good semiquantitative agreement. For both experiment and theory normal state spectra of nearly 
optimally doped PCCO show clear evidence for a pseudogap state with AFM-like nature. Namely, 
folding of quasiparticle bands as well as presence of the "hot spots" and "Fermi arcs" were observed. 

PACS numbers: 74.72.-h; 74.20.-z; 74.25. Jb; 31.15.A- 



I. INTRODUCTION 

Many experimental and theoretical papers have been 
seeking to describe the nature of high-temperature super- 
conductivity (HTSC) in cuprates. In contrast to the nor- 
mal (Fermi-liquid) metal, HTSC compounds show many 
abnormal properties for temperatures above the super- 
conducting transition; the normal state pseudogap being 
a notorious example 1 . The origin of this anomalous state 
is usually attributed either to superconducting fluctua- 
tions (precursor Cooper pairing) 2 or to some order pa- 
rameter competing with superconductivity 3,4 , e.g. AFM 
fluctuations, incommensurate or fluctuating charge den- 
sity waves (CDW), stripes, etc. 

Recently a generalized LDA+DMFT+Sk computa- 
tional scheme was proposed to describe the pseudogap 
state in strongly correlated systems, by accounting for 
nonlocal AFM (or CDW) fluctuations with short-range 
order 5 ' 6 ' 7 . Its relation to other theoretical DMFT-likc 8 



approaches to the pseudogap state was discussed e.g. in 
Ref. 4. Both model computations and those for real 
systems were done 4 ' 9 . This approach, for instance, al- 
lowed one to describe the experimentally observed par- 
tial Fermi surface (FS) "destruction" 5 , which was the- 
oretically studied for hole-doped HTSC prototype sys- 
tem Bi 2 Sr2CaCu208-(5 (Bi2212) 4 and electron-doped one 
Ndi.85Cerj.i5Cu04 (NCCO) 9 . Two-particle properties 
can also be described by this approach 7 , e.g. calculated 
optical spectra in the pseudogap state compare well with 
experimental data for Bi2212 4 and NCCO 9 . 

In this paper we study the electron-doped 
Pr2-a;Cc K Cu04 (PCCO) in the pseudogap state 
(x = 0.15) using the generalized LDA+DMFT+S k com- 
putational scheme 5 ' 6 ' 7 and ARPES measurements 10,11 . 
We present here both experimental and theoretical 
quasiparticle spectral functions and Fermi surfaces. 
These are found to agree well with each other supporting 
competing order parameter fluctuations as origin of the 
pseudogap instead of superconducting scenario. 
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II. COMPUTATIONAL DETAILS 

The crystal structure 12 of P^CuCU has tetragonal 
symmetry and the space group is I4/mmm. Lattice 
constants 12 are a=b=3.962, c=12.154 A. There are two 
crystallographic types of oxygen atoms in P^CuO^ the 
first one belongs to CuC>2 layer and the second is located 
within the Pr atoms. The atomic positions in the ele- 
mentary cell are as follows: Cu - 2a(0,0,0), 01 
- 4c(0,0.5,0), Pr - 4e(0,0,0.35171), 02 - 4d(0,0.5,0.25). 12 




FIG. 1: The crystal structure of Pr2Cu04. A middle size 
grey spheres correspond to the copper atoms, a small dark 
and black spheres are Ol and 02 atoms, respectively and a 
big grey spheres — praseodymium atoms. 

In Fig. 1 the crystal structure of Pr2Cu04 is shown. 
Middle size grey spheres correspond to the copper atoms, 
little dark and black spheres represent 01 and 02 atoms, 
big grey spheres show praseodymium positions. Clearly 
visible quasi two-dimensional nature of these compounds 
determines its physical properties. Physically most in- 
teresting are the Cu02 layers. Those layers provide an- 
tibonding Cu-3d(a; 2 — y 2 ) partially filled orbital, whose 
dispersion crosses the Fermi level. Thus we are using this 
effective Cu-3c?(x 2 — y 2 ) antibonding band as a "bare" 
band in DMFT computations. 



> 



r 

(0:0) 



x 

too) 



M 

Cnpr) 



r 

(0;0) 



FIG. 2: LDA+DMFT+S fc quasiparticle energy dispersion of 
PCCO Cu-3d(a; 2 — y 2 ) orbital for high symmetry directions 
of square Brillouin zone. The Fermi level is zero. 



Electronic structure of PCCO was investigated within 
generalized LDA+DMFT+E k scheme. 5 ' 6 ' 7 This scheme 
incorporates the density functional theory in local den- 
sity approximation (LDA) 13 and the dynamic mean- field 
theory (DMFT) 8 with "external" momentum-dependent 
self-energy Ek 6 . 

As a first step the LDA band structure calculation 
was performed. Using crystal structure data, the elec- 
tronic band structure was obtained with the linearized 
muffin-tin orbitals (LMTO) method 14 . Further hopping 
integral values were calculated for effective Cu-3<i(a; 2 — 
y 2 ) Wannier function within the N-th order LMTO 
(NMTO) framework. 15 Corresponding hopping integral 
values are t = -0.4385, t' = 0.1562, t" = 0.0976. The 
value of Coulomb interaction on effective Cu-3g?(x 2 — y 2 ) 
orbital 17=1.1 eV was obtained via constrained LDA 
computations 16 . Secondly, the DMFT calculations 8 , 
which take the hopping integrals and the Coulomb in- 
teraction as input parameters, were performed. 

To account for the AFM spin fluctuations, a two- 
dimensional model of the pseudogap state is applied. 17 
Corresponding k-dependent self-energy Ek 1 ' 17 describes 
nonlocal correlations induced by (quasi) static short- 
range collective Heisenberg-like AFM spin fluctuations. 
The quasi static approximation for AFM fluctuations 
necessarily limits our approach to high - enough tem- 
peratures (energies not so close to the Fermi level) 17 , so 
that in fact we are unable to judge e.g. on the nature of 
low temperature (energy) damping in our model. Thus 
we avoid here possible discussion of whether the damping 
corresponds to marginal or regular Fermi liquid. More- 
over as shown in Ref. 6 for the "hot-spot" corresponding 
self-energy has essentially non Fermi liquid behavior. 

The Ek definition contains two important parameters: 
the pseudogap energy scale (amplitude) A, representing 
the energy scale of fluctuating SDW, and the spatial cor- 
relation length £. The latter is usually determined from 
experiment. The A value was calculated as described in 
Ref. 6 and found to be 0.275 eV. The value of correlation 
length was taken to be 50 lattice constants, in accordance 
with the typical value obtained in neutron scattering ex- 
periments on NCCO 18 . To solve DMFT equations nu- 
merical renormalization group (NRG 19,20 ) was employed 
as an "impurity solver" . Corresponding temperature of 
DMFT(NRG) computations was 0.011 eV and electron 
concentration was n=1.145. 



III. EXPERIMENTAL DETAILS 

Photocmission experiments were performed at UE112- 
PGM beamline at BESSY using SCIENTA SES100 
analyzer. Typical energy and angular resolution for 
the excitation energy (hv = 100 eV) used in this 
study were 20 meV and 0.2° respectively. Samples of 
P r i.85Ceo.i5Cu04 +( 5 were grown using traveling solvent 
floating zone technique and annealed to achieve optimal 
T c of 25K with transition width of IK, which resulted 
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FIG. 3: (a) Extended Fermi surfaces for PCCO - 
LDA+DMFT+Ek data. White rectangle on panel (a) 
schematically shows the part of reciprocal space measured 
experimentally (panel b). Lower left corner is X-point (tv, 0). 



from improved growth conditions 21 . Similarly the with 
(FWHM) of X-ray rocking curves was less than 0.08°, 
signaling high quality of the samples. 

For the photocmission measurements the samples were 
mounted on a cryomanipulator and cleaved in situ in ul- 
tra high vacuum with a base pressure < 1 • 10~ 10 mBar. 
During the whole experiment, including the temperature 
cycling, when the sample was heated to room temper- 
ature and then cooled back again, no observable aging 
effects were detected. 



folding, as it would be the case for a long-range AFM or- 
der. Namely, quasiparticle band and the "shadow" band 
are not exactly the same. No real band gap opens at 
(tt/2, 7r/2) point. Nevertheless, suppression of the spec- 
tral weight is clearly detectable in the vicinity of X (ir, 0) 
point, thus signaling opening of the pscudogap, which in 
this case can be viewed as a precursor of the real band 
gap. Splitting takes place between the quasiparticle band 
and the "shadow" band with the value of about 2A. 

In Fig. 3 an extended picture of PCCO Fermi sur- 
faces is presented (panel (a) — LDA+DMFT+Ek results, 
panel (b) — experimental ARPES data). Strictly speak- 
ing Fig. 3 is a color map in reciprocal space of the cor- 
responding spectral function plotted at the Fermi level. 
FS is clearly visible as reminiscence of non-interacting 
band close to the first Brillouin zone border and around 
(tt/2, 7t/2) point (so called Fermi arc), where the quasi- 
particle band crosses the Fermi level. There is an inter- 
esting physical effect of partial "destruction" of the FS 
observed in the "hot spots", points that are located at 
the intersection of the FS and its AFM umklapp replica. 
This FS "destruction" occurs because of the strong elec- 
tron scattering on the antifcrromagnctic (AFM) spin 
(pseudogap) fluctuations on the copper atoms. Also the 
"shadow" FS is visible as it should be for AFM fold- 
ing. As no long-range order is present in the underdoped 
phase the "shadow" FS has weaker intensity with respect 
to FS. Another evidence of presence of electron pock- 






IV. RESULTS AND DISCUSSION 

Generally speaking, finite temperature and interaction 
lead to notable life-time effects. Thus, instead of quasi- 
particle dispersions obtained in usual band structure cal- 
culations one has to deal with corresponding spectral 
function A(ui, k): 

A(u,k) = --lmG(u,k), (1) 

where G(w,k) is the retarded Green's function obtained 
via LDA+DMFT+Ek scheme 5 ' 6 ' 7 . 

Color plots, whose intensity encodes the function val- 
ues became a traditional and convenient way of represent- 
ing these multiple variable functions. Such a color plot of 
LDA+DMFT+Ek quasiparticle spectral function (1) of 
copper 3d(x 2 — y 2 ) orbital is presented in Fig. 2. Width 
of the quasiparticle spectral function in the color plot is 
inversely proportional to the quasiparticle life time. The 
calculated quasiparticle band dispersion has minimum at 
r point (-1.52 eV) and maximum at M point (2 eV). It 
crosses the Fermi level along X-M as well as M-r direc- 
tions. Because of AFM pscudogap fluctuations there is a 
well detectable (quasi) folding of the quasiparticle band, 
reflected in the formation of the so called "shadow" band, 
which has its maxima at T-point and minima at M-point. 
However, because of the short-range nature of the anti- 
ferromagnetic order, this does not result in a complete 
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FIG. 4: Energy-momentum intensity distributions for the 
specific cuts drawn in Fig. 3 (upper panels — theoretical data, 
lower panels — experimental photoemission intensity). To 
judge about the absolute intensities of the "shadow" (1) and 
main band (2) cut 1 contains an MDC curve integrated in 
an energy window 60 meV centered at the Fermi level (FL). 
Similarly integral EDC for cut 2 ( "hot spot" ) shows suppres- 
sion of the intensity at the FL as compared to cut 3, which is 
located further away from the "hot spot" . The FL is zero. 
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ets can be seen in the experimental FS map shown in 
Fig 3b. One pocket is centered at point with coordinates 
(0,0. 8ir) and the other at (0.87r,0). Again, in agreement 
with the theoretical prediction, the pocket "sides" which 
coincide with originally unreconstructed FS have higher 
intensity, while the newly appearing "replicas/shadows" 
have weaker intensity. 

The PCCO FS is very similar to that of 
Nd2- 2; Ce a ;Cu04 (NCCO), which belongs to the same 
family of suppcrconductors. The NCCO was recently 
studied both theoretically 9 and experimentally 22 . 

Let us compare theoretical (upper line) and experimen- 
tal (lower line) energy quasiparticle dispersion for most 
characteristic cuts introduced in Fig. 3 (see Fig. 4). The- 
oretical data were multiplied by the Fermi function at 
a temperature of 30K and convoluted with a Gaussian 
to simulate the effects of experimental resolution, with 
further artificial noise added. 

The Cut 1 intersects quasiparticle and "shadow" Fermi 
surfaces close to the Brillouin zone border. One can 
find here a "fork" -like structure formed by the damped 
"shadow" band (-0.5-0 arb. u.) and better defined quasi- 
particle band (0.5-1 a.u.). This structure corresponds to 
preformation of FS cylinder around (7r,0) point. The Cut 
2 goes exactly through the "hot spot". Here we see a 
strong suppression of the quasiparticle band around the 
Fermi level as it is also shown in Fig. 2. The Cut 3 crosses 
the Fermi arc, where we can see a very well defined quasi- 
particle band. However weak intensity "shadow" band 
is also present. For the case of long range AFM order 
and complete folding of electronic structure, FS and its 
"shadow" should form a closed FS sheet around (ir/2, 
it/2) point, while in the current case the part of the 
pocket formed by the "shadow" band is not so well de- 
fined in momentum space. As can be seen there is a good 



correspondence between the calculated and experimental 
data in terms of the above described behavior, which is 
also similar to the results reported for Nd2- x Ce x Cu04 
(NCCO) in our earlier work. 9 



V. CONCLUSION 

In this work the LDA+DMFT+S k was performed for 
electron-doped Pr2_2;Ce x Cu04 compound in the pseu- 
dogap regime. The LDA+DMFT+Ek calculation shows 
that Fermi-liquid behavior is still conserved far away 
from the "hot-spots" , while the destruction of the Fermi 
surface observed in the vicinity of "hot spots" is due to 
strong scattering of correlated electrons on short-range 
antifcrromagnetic (pseudogap) fluctuations. Comparison 
between experimental ARPES and LDA+DMFT+E k 
data reveals a good semiquantitative agreement. The 
experimental and theoretical results obtained once again 
support the AFM scenario of pseudogap formation not 
only in hole doped HTSC systems 4 but also in electron 
doped ones 9 . 
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